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Abstract 

We use Lorentz polynomials to present the solutions explicitly of equations 
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1 Introduction 

It is well-known that Multiresolution Analysis (MRA) is a general procedure to 
construct wavelet basis. 

Definition 1.1 A Multiresolution Analysis (MRA) generated by function (j) con¬ 
sists of a sequence of closed subspaees Vj,j € h, o/L 2 (R) satisfying 

(i) (nested) Vj C g+i for all j G Z; 

(ii) (density) UjgzP = r 2 (R); 

(iii) (separation) fljgzVj = {0}; 

(iv) (scaling) f{x) G Vj if and only if f{2x) G V)+i for all j G Z; 

(v) (Basis) There exists a function (f G Vq such that {4>{x — k) : k G K) is an 

orthonormal basis or a Riesz basis for Vq . 
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The function whose existence asserted in (v) is called a scaling function of the MRA. 

A scaling function (f must be a function in ^ 2 ( 1 ^) with f if ^ 0. Also, since cf G Vq 
is also in Vi and {(fi^k ■= 2^/^(f>{2x — fc) : fc G Z} is a Riesz basis of Vi, there exists a 
unique sequence {pfe}^_oo G hi'^) that describe the two-scale relation of the scaling 
function 


X! Vk<t){2.x-k), (1.1) 

k— — oo 

i.e., is of a two-scale refinable property. By taking a Fourier transformation 
on both sides of (1.1) and denoting the Fourier transformation oi f by := 

(l>{x)e~'^^^dx, we have 


where 


( 1 . 2 ) 


1 

P{z) = - ^ pkz’^ and z = (1.3) 

k— — oo 

Here, P{z) is called the mask of the scaling function. Now, regarding the property 
that {(j){x — k)} must be an orthonormal basis, we have the following characterization 
theorem (see, for example, Chs. 2, 5 and 7 of Chui [5] and Ch. 3 of Hernandes and 
Weiss [7]) 


Theorem 1.2 Suppose the function (f satisfies the refinement relation (fix) = 
~ k)- Then (i) {(f{x — k) : k G mathbbZ} forms an orthonormal basis 
only if \P{z)\‘^ + \P{—z)\'^ = 1 for z G C with \z\ = 1. 

(a) Suppose P{z) satisfies 

1. P{z) G and is 2Tr-periodic 

2 . |P(z)|2+|P(-z)|2 = l 

3. P(l) = 1 

4 . P{z) ^ 0 for all ^ G [—tt, tt] 

Then {if{x — k) : k G mathbbZ} forms an orthonormal basis. 


Denote the cardinal B-splines with integer knots in No by Bn{x). It is well-known 
that Bn{x) satisfy rehnement relation 
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S„(x) — ^ 


i=o 


Bn(2x - j) 


(1.4) 


and have masks Pn{z) such that Bn{0 = Pn{z)Bn{i/2), where z = e and 


1 ” 1 

Pn{z) = 2 X! 2^ 
j=o 


{l + zY 


1 + z 


(1.5) 


It is clear that 

\Pn{zr + \Pn{-z)f = 


1 + Z 


2n 


1-Z 


2n 


= cos^"(^/4) + sin^"(^/4) < cos^(^/4) + sin^(^/4) = 1 


The equality happens only when n=l. Therefore, except for the case of order one 
(i.e., n = 1), Bn{x) are generally not orthogonal (indeed they are Riesz basis). To 
induce orthogonality, Daubechies (see [11 [5]) introduces a class of polynomial function 
factors S{z). Hence, instead of Bn{x), a scaling function (j)n{x), called spline type 
scaling functions, with the mask Pn{z)Sn{z) is considered so that 


(j)n(.C) = Pn{z)Sn{z)(j)n{^/2), ( 1 . 6 ) 

where P„(z) are defined as (ILSI). We need to construct S'„(z) such that the shift set 
of the new scaling function form an orthogonal basis. In other words, we need that 
Sn{z) satisfy the following condition 

\Pn{z)Sn{z)\^ + |P„(-z)5„(-z)P = 1 (1.7) 

Now we consider Sn{z) of the following type: Sn{z) = 012 : + 022 ^ + ... + a„z”, n G N 
and Oi G M, i = l..n. When 2 ; = 1, from equation (II3 we have 

1 = |Pn(l)5n(l)|" + |Pn(-l)5„(-l)P 
= ip„(i)n5„(i)p + ip„(-i)n5„(-i)p 
= |5„(l)p + 0 = |5„(l)p 

Thus S'„(l) = X)r=i From Theorem 1.2 (ii), we further impose a re¬ 

striction that X)i=i Oi = 1 in order to ensure the orthogonality of the scaling function. 

Next, we set out to find the expressions and constructions of Sn- We have the 
following Lemma and leave the proof for next section. 
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Lemma 1.3 

Let Sn{z) be defined as above. Then there holds 


\Sn{zr 

n n—1 n—2 

= X! ^ X! a^ai+lCos{£,/2) + 2 ^ aia^+2Cos{2^/2) + ... 

i—1 i—1 i—1 

+ 2aianCos{{n — 1)^/2). 


From Lemma 1.3, if we write each cos{k^/2) as a polynomial of cos(^/2), then 
\Sn{z)\^ = Qn{x) where x = cos{£^/2). Obviously, Qn{x) has the degree of n — 1. It 
is also easy to observe that |S'„(—z)p = Qn{—x). Now equation (2.5) becomes 


1 = \Pn{z)Sn{z)\^ + \Pn{-z)Sn{-z)\^ 

= cos^"(5/4)Q„(a:) + sin^”(^/4)Q„(-a;) 


1 + cos{£,/2) 


1 + X 


Qn{x) “b 


Qn{x^ “b 

1 — X 


1 — cos{^/2) 

2 

n 

Qni^ 


Qn(^ 


So finally we get 


(^-^) Qn{x) + '5n(-2;) = 1 (1-8) 

As a side note, (11.81) is equivalent to (6.1.7) of [4] or (4.9) of [5], but is of quite 
different form so that we may obtain a complete different solution of the equation 
shown below in (HU by using Lorentz polynomials, which yields an efficient proof of 
sufficiency for the orthogonality and a mechanical and elementary way to construct 
scaling functions (fn- 


Next, to show the existence of Q{x) in the above equation, we make use of the 
Polynomial Extended Euclidean algorithm (see Cormen, Leiserson, Rivest, and Stein 

[3])- 


Lemma 1.4 Polynomial extended Euclidean algorithm If a and b are two 

nonzero polynomials, then the extended Euelidean algorithm produces the unique pair 
of polynomials (s, t) such that as+bt=ged(a,b), where deg{s) < deg(h)—deg(gcd(a, b)) 
and deg(t) < deg{a) — deg{gcd{a,b)). 

We notice that (?cfi((i^)”, (^^^)") = 1, so by Lemma 1.4, there exists uniquely Q{x) 
and R{x) with degrees less than n such that (^^)^Q{x) + {^-^)^R{x) = 1. If we 
replace x by —x in the previous equation, we have (^!^)"(5(—x) + (^!^)”i?(—x) = 1. 
Due to the uniqueness of the algorithm, we conclude that i?(x) = Q{—x). So we 
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have showed the existence of a unique Q(x) = Qn{x) satisfying equation (11.81) . 


To construct Qn(x) explicitly, we use the Lorentz polynomials shown in Erdelyi 
and Szabados [6], Lorentz [10] and the following technique. 



where the polynomials presenting in the brackets are the Lorentz polynomials. We 
notice that the degrees of the two polynomials in the brackets are n — 1, and because 
Qn{x) in equation 2.6 is unique, we can conclude that 



With the construction of Qnix), we take a step further by showing the existence of 
S S fliO-i+ii ■•■in Lemma 1^3. 

It is well-known that the set {l,cos{t),cos{2t), ...,cos{(n — l)t)} is linearly inde¬ 
pendent. As a result, {1, cos(^/2), cos(2^/2),..., cos((n— 1)^/2)} forms a basis of the 
space Pn-i{x) = {P{x) : x = cos{^/2) and P is a polynomial of degree less than n. 
Based on this fact and the existence of Qnix) in equation (2.6), it is obvious that 
the coefficients •■■ in Lemma 2.2 must exist uniquely. 

We now survey the results of $6.1 of |1] and $4.B. of |S] in the following theorem 
and leave its proof by using our mask expression (HH) for next section. 

Theorem 1.5 Let P„(z) and Sniz) be defined as above. Then for n = 1, 2,..., 
the spline type function finix) with the mask Piz) = Pniz)Sniz), where Sniz) = 
aiz + a 2 z‘^ -b ... -b UnZ^, is a scaling function that generates an orthogonal basis of Vq 
in its MR A. 

Although we can not give the explicit expression of finix), we may use the recur¬ 
sive method presented in Theorem 5.23 of [1] by Boggess and Narcowich to find 
an approximation of finix) with any accuracy. It is easy to see that all three 
conditions required by Theorem 5.23 are satisfied by Pniz): (i) Pn(l) = 1, (h) 
|Pn(-z)P + \Pni—z)\‘^ = 1 (|.z| = 1), and (iii) Pniz)\ > 0 (5m[— tt,tt]). The examples 
for the cases n = 3 and 4 will be provided in Section 3 to demonstrate this procedure, 
while all proofs are given in next section. 
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2 Proofs 

Proof of Lemma 1.3 We have 
\Sn{z)\^ 

= |ai(cos(^/2) — isin{^/2)) + ... + an{cos{n£^/2) — isin{nf^/2))\‘^ 

= |(aicos(^/2) + ... + a„cos(n^/2)) — i{aisin{^/2) + ... + a„szn(n^/2))p 
= (aicos(^/2) + ... + a„cos(n^/2))^ + (aiSzn(^/2) + ... + a„szn(n^/2))^ 
= ai(cos^(^/2) + szn^(^/2)) + ... + a„(cos^(n^/2) + sm^(n^/2)) 

+ 2aiaj{cos{if,/2)cos{j^/2) + sin{i^/2)sin{jf./2)) 


n 

= ^ ^ 2aiajCos{{i — j)^/2) 

1=1 iyj 

n n—1 n—2 

= X! ^ X! ajai+icos(C/2) + 2 ^ ajaj+2Cos(2^/2) + ... 

i—1 i—1 i—1 

+ 2aia„cos((n — 1)^/2) 

A similar procedure can be applied to find |S'„(—z)p 

|54-^)r 

n n—1 n—2 

a^ai+icos{^/2) + 2 ^ a^a^+2Cos(2^/2) + ... 

i—1 i—1 i—1 

+ (-l)”2aia„cos((n - l)C/2). 


We now prove Theorem 1.5 using our new mask expression (HU). 

Proof. From theorem 1.2 (ii), the sufficient conditions for the orthogonality of the 
scaling function are 

1. P{z) G and is 27r-periodic 

2. |P(z)|2+|p(_^)|2 = l 


3. P(l) = 1 

4. P{z) ^ 0 for all ^ G [—tt, tt] 


From the construction of our P{z) = Pn{z)Sn{z), the first two conditions are au¬ 
tomatically satisfied. The third condition is also obvious: P(l) = Pri(l)<5'„(l) = 
-) X]r=i = 1 according to the construction of Sn{z). Now we will prove that 


' 1+1 ^ 


the final condition is fulfilled as well. 
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Indeed, if C S [—7r,7r], then firstly we have 
\Pniz)\ = 


1 + 

n 

1 + cos{f/2) — isin{f/2) 

2 


2 


= |cos^(^/4) — ism(^/4)cos(^/4)|" = \/cos^^f7^T+'cos^(^7^05^^^^(^74) 
= |cos(^/4)|" > |cos(7r/4)|" > 0 for ^ G [—tt, tt] 

Secondly, from equation (HH) we have 



Thus IS'nC-z)! > 1 and \P{z)\ = \Pn{z)\\Sn{z)\ > |cos(7r/4)|” > 0. 

The remaining thing is to show that the newly constructed scaling function (x) 
with mask Pn{z)Sn{z) is in 4^2(2), where Sn{z) = n G N. From He 

019], we know that G L 2 (M) if 

n 

( 2 . 10 ) 

i=i 

Recall from Lemma 1.3 that 

Qn{x) = QnicOs{f,/2)) 

n n—1 

= E«' + 2E aiai+icos(^/2) + ... + 2aia„cos((n — l)C/2) 

2—1 2—1 

Taking the integration from 0 to 27r of both sides, we have 



Qn{x')d^ 



aiai+icos{f,/2) + ... + 2aia„cos((n — 



= 27r^a? 


df 
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On the other hands, from the expression of Qn{x) in we have 




Qn{x)d^ = 


2 = 0 




f 'e (i d( 


Combining the two equations above we have 


10 -j p27T 1 

E“? = ^/ E 


2=1 


2=0 


2n — 1 \ /1 + cos(^/2) 


n—1 —2 


1 — cos{^/2) 


d^ 

( 2 . 11 ) 

Now we need to show that the expression on the right hand side of (12.1111 is smaller 
— 1 

than - -. 


First of all, it is easy to see that for 0 < i < n — 1 


2n - 1 \ / 2n - 1 \ _ 1 / 2n 

i /^In — 1/ 2vri 


Applying this inequality into (12.lip yields 


J_ 2n - 1 

27r ^ 


i=0 


/O 


l + cos(^/2)\" ^ */l —cos(^/2) 


di 


i=0 "'o 


47r V « / ./n V 4 

^ 2 — r 

"-1 fTr/2 


1 / 2n 


TT \ n 


^ ^TT/Z 

(sma;)^*dx 

2=0 -^0 


for a: = ^/4. 
Now let 


n —1 


—^ /•7r/2 

2=0 -^0 

We can express A as the sumation of two terms A = Ai + A 2 , where 
Ai = 'Y / {cosx)'^^^~^~^\sinx)'^^dx 

.■_n Jo 


2=0 


n-l pTT/2 

A2 = Y^ / (cosa;)^^"“^“*ysinx)^*(ia; 

i=[2^i] + l ° 


( 2 . 12 ) 
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For 0 < i < 2(n — 1 — i) > 2i, the term Ai becomes 


1 —l l 

2 J /*'n'/2 


^ / {cosx)'^^"'~^~'^^\sinxcosx)'^^dx 

.■_n Jo 


= ^ — J (cosa;)^^" ^ ‘^^\sin 2 x)^''dx 




For + l<z<n — 1, 2(n — 1 — z) < 2z, the term ^2 becomes 

n-l „^/2 


^2 = ^ / (sznxcosx)^^"“^“*^(sznx)^^^*“"^^^dx 

i=[^] + l ° 

"-1 1 ^7r/2 


i=[2^] + l 


^ E 


I C’r/2 


(szrzx)"(2*-"+i)rf^ 




Next, we make use of the following well-known result 


■ ^ 2 n. r'\ ^ 2 r^. (2n-l)!!^ 

/ (stnx) dx= (cosx) dx = —;——- 

Jo Jo (2n)!! 2 


(2n)! TT TT (2n)! tt 1 / 2n 
[(2n)!!]2 2^2 4"(n!)2 ^ 24" V ^ 


Next, we try to find the upper bound for 
[16], we have the following inequalities 


. Based on Stirling estimation in 


/3n+T 


Tt } x/TTfi V 12n — 1 
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Using (3.6) on Ai and A 2 yields 


I 2 I 1 

^ = ^1+^2 <2 


1—1 , 
=^y I- 


^ 4* ^3(n - 1 - 2z) + 1 2 


1 


^ 3* (|)V3(ri-l-2*) + l 
We consider the denominator of the fraction, and let 

A on _ 1 

fix) = i^)^^[Hn -l-2x) + l],xe [0, 

Surveying the function, we have f(x) attains minimum at 0, or 

fix) = (^)'"[3(n - 1 - 2a;) + 1] > /(O) = 3(n - 1) + 1 

Thus, we have 


A<n^ 


1 1 


< 


00 

TT 1 


y- = 

Ql 


Stt 


3 * y'3(n-l) + l “ \/3^r^ ^ 3* 2y/3^fr^ 

Finally, combining (I2.11I) - (I2.15|1 . we have 

1 pn 1^/ 1__^\ 3^ 

^ TT \ n J tt y/im \ 12n — 1 / 2-v/3n — 2 


(2.15) 


<^2n — 1 

For n > 17, we can easily verify that the right hand side is less than - -. Using 

Mathematica for direct calculation of the case n < 16, we Hnd that the inequality 
in theorem 3.1 holds. Thus, it holds every interger n, and we have shown that 
(j) S L 2 (R) and complete the proof. 


3 Examples 


It is easy to find that 51 ( 2 ;) = 2 and (j)nix) is the Haar function. For n = 2, 


5 * 2 ( 2 ) 


1 + ^3 l-\/3 2 

—— 2 ; H--— 2 "^ 


and the corresponding (f^ix) is the Daubechies D 2 scaling function. 

As examples, we consider the efficiency of the computation in using mask ex¬ 
pression (fm to construct Daubechies scaling functions (the spline type scaling 
functions) for the cases of n = 3 and 4. According to the previous sections, we will 
construct the Daubechies scaling function (/> 3 (x) from the third order B-spline 
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function B'i{x) using our expression (11.91) . 


In order to construct the function (j)^{x), we start with its mask P 3 (z)S 3 {z), where 
P 3 {z) = is the mask of the third order B-spline. Let S 3 {z) = aiz+a 2 z‘^+a 3 z^, 

then by Lemma 1.3, we have 

Qaix) = |5'3 (z)P = + al+ Og) + 2(0102 + 0203)005(^/2) + 2oi03Cos(^) 

= (o^ “t“ O2 Og) 2 (oi 02 “t“ 020g)c05(^/2) — 1) 

= (oi + of + Og - 201O3 ) + 2(oi02 + 0203)005(^/2 ) + 4oi03COS^(^/2) 

= (of of “f O 3 — 2 O 1 O 3 ) p 2 (oi 02 “t“ 0203 )^ p 4 oi 03 X^ 

where x = cos(^/2) and z = 


On the other hand, by equation dm) we have 


< 33 ( 2 ;) = ^ 


i =0 


IP X 


\P X 


3 2 9 

= — -X 


+ 5 
+ 4 


1 — X 
2 

1 — X 


+ 10 


1 — X 


2 2 

Thus, from the above equations, we have the following system of equations 


( of + of + of — 201 O 3 = 4 

2(0i02 + O2O3) = -| 


(3.16) 


[4oi03 = 2 


Simplify (4.1) we get 


of + of = f 
01 O 2 + 02 O 3 = —I 


(3.17) 


0103 — g 

From this system, we have (oi + 02 + 03)^ = of + of + of + 2 (oi 02 + 02O3 + 01O3) = 1. 
Without loss of generality, consider the case oi + 02 + 03 = 1. Combining this with 
01 O 2 + 02 O 3 = —I and 01 O 3 = I we have the following solution 


(3.18) 


' = 3 (1 + \/T0- \/5 + 2^) 

02 = i (1 - yio) 
as = 3 (1 + yiO + v'S + 2^) 


We verify the condition for 1/3 to be in L 2 { 


2 2 19 32 22-3-1 

^2 + «3 — ^ —- 


3 
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Thus (f>3 is indeed in L2(K.). Now we will attempt to construct (f>3 explicitly. It has 
the mask 

P3{z)S3{z) = ^ (01^ + 022:^ + 032;^) 

= 0.0249a; - O.OOOdx^ - O.OOSai^ + 0.325a;^ + 0.571x® + 0.2352a;® 
Hence, we have the refinement equation 

((13(2;) = O.O 49803 ( 2 a; - 1 ) - 0.121((i3(2a; - 2 ) - 0 . 191 (/> 3 ( 2 a; - 3 ) 

+ 0 . 650 ((i 3 ( 2 a;- 4 ) + 1 . 141 ^ 3 ( 2 a;- 5 ) + 0 . 4705 ^ 3 ( 2 a;- 6 ) ^ ' 

Again, we construct the Daubechies D4 orthogonal scaling function (j)4{x) from 
the forth order B-spine i?4(x) with the mask ^4(2;) = ■ 

Now we examine the mask P4(z)5'4(z) of (j)4{x) where we define S'4(2:) = 012: + 
022:^ + 03Z® + 042:"'’. By Lemma 1 . 3 , we have 

<34(2;) = |<S'4(2:)P = (oi + 02 + 03 + 04) + 2 (oi02 + 02O3 + 0304)005(^/2) 

+ 2(0103 + 0204 )cos(^) + 2oia4Cos(3^/2) 

= (oi + O2 + O3 + O4) + 2 (oi02 + 02O3 + 0304)005(^/2) 

+ 2(0103 + a204)(2oo5^(^/2) — 1 ) + 2 oi 04(4005^(^/2) — 3 oo 5 (^/ 2 )) 

“ (Oi + O 2 + O 3 + O 4 — 2 oi 03 — 20204 ) 

+ (2 oi 02 + 20203 + 20304 — 60104)005(^/2) 

+ (4oi03 + 40204 ) 005 ^(^/ 2 ) + 80104005 ^(^/ 2 ) 

= (Oi + O2 + O 3 + O 4 — 201O3 — 202O4) 

+ (2oi02 + 202O3 + 203O4 — 60104)0; 

+ (4oi03 + 40204)2;^ + 801043;® 

where x = cos[^/2) and 2; = 


Using equation (EH we get another expression for < 34 ( 2 ;) 



Thus, from the above equations, we have the following system of equations 

of + O 2 + a| + a| — 20103 — 20204 = 8 
2 ai 02 + 20203 + 20304 — 60104 = — ^ 

4ai03 + 40204 = 10 
80104 = — I 


(3.20) 
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Simplify (5.1) we have the following system 


af + a| + 03 + 04 = 13 

1 SI 

0102 + 0203 + O3O4 = — ^ 
O1O3 + O2O4 = I 

5 

CL1CL4 — yF 


(3.21) 


Solving for this system of equations yields 8 solutions. One of the numerical solutions 
is 

( 01 = 2.6064 


02 = -2.3381 
03 = 0.8516 
04 = -0.1199 


(3.22) 


We verify the condition for to be in L 2 (K) 


Oi -f O 2 “t“ O 3 -\- O 4 — 13 < 32 — 


22.4-1 

4 


Thus 04 is indeed in L 2 (®.). Now we will attempt to construct 04 explicitly. It has 
the mask 


Pi{z)SA{z) = + asz^ + 04 z"‘) 

= 0.1629Z + 0.5055^2 + 0.44612^ - 0.0198z'‘ - 0.1323z^ + 0.0218z® 
+ 0.0233z’^ - 0.0075Z® 


Hence, we have the refinement equation 

04(a;) = O.325804(2 x - 1) + l.O1104(2a; - 2) + O.892204(2a; - 3) - O.O39604(2x - 4) 

- 0 . 264604(21 - 5) + O.O 43604 ( 2 a; - 6 ) + O.O46604(2a; - 7) - O.O 1504 ( 2 a; - 8 ) 

(3.23) 

The above examples demonstrate the efficiency of the computation of the Daubechies 
scaling functions by using expression (fm . 
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